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ABSTRACT 

Middle-aged, cooling neutron stars are observed both as relatively rapidly spinning radio pulsars 
and as more slowly spinning, strongly magnetized isolated neutron stars (INSs), which stand out by 
their thermal X-ray spectra. The difference between the two classes may be that the INSs initially had 
much stronger magnetic fields, which decayed. To test this, we used the Chandra X-ray Observatory 
to observe IRXS J072559.8-261229, a possible X-ray counterpart to PSR J0726-2612, which, with 
its 3.44 s period and 3 x 10^'^G inferred magnetic field strength, is the nearest and least cxtincted 
among the possible slowly-spinning, strong-field INS progenitors (it likely is in the Gould Belt, at 
~ 1 kpc) . We confirm the identification and find that the pulsar has a spectrum consistent with being 
purely thermal, with blackbody temperature kT = 87±5eV and radius R = 5.7^^'3km at a distance 
of 1 kpc. We detect sinusoidal pulsations at twice the radio period with a semi- amplitude of 27 ± 5%. 
The properties of PSR J0726— 2612 strongly resemble those of the INSs, except for its much shorter 
characteristic age of 200 kyr (instead of several Myr). We conclude that PSR J0726— 2612 is indeed 
an example of a young INS, one that started with a magnetic field strength on the low end of those 
inferred for the INSs, and that, therefore, decayed by a relatively small amount. Our results suggest 
that the long-period, strong-field pulsars and the INSs are members of the same class, and open up 
new opportunities to understand the puzzling X-ray and optical emission of the INSs through radio 
observations of PSR J0726-2612. 

Subject headings: stars: individual (PSR J0726— 2612) — stars: neutron — X-rays: stars — X-rays: 
individual (IRXS J072559.8-261229) 



1. INTRODUCTION 



The ROSAT All-Sky Survey fRASS : IVoiii"ird1[l999l) 

showed that our census of cooling, nearby neutron stars 
was incomplete: it contained not just the known cool- 
ing pulsars such as PSR B0656+14 and Geminga, but 
also seven "isolated neutron stars" (INSs). These are 
nearby (< 1 kpc), young (< 1 Myr), cooling neutron stars 
with thermal X-ray spectra, long periods (> 3 s), faint 
optical coun terparts, and no detected radio emission (for 
reviews, see iHaberll [20071 iKaplanI [20081 ). The INSs are 
interesting both because of their abundance and because 
of the promise of inferring neutron-star parameters from 
their thermal emission. Unfortunately, despite large in- 
vestments of time with Chandra and XMM-Newton, the 
nature of the emission remains puzzling, and we still un- 
derstand neither the composition nor state (gaseous, con- 
densed) of the surface. 

A clearer picture has emerged for the origin and abun- 
dance of the INSs: most likely, they were born with 
very strong magnetic fields, of > lO^"^'^ G, which de- 
cayed. Empirical ev idence for this comes from o ur 
X-ray ti ming e fforts (Kaplan & van Kerkwiiy[2005iB, 
I2009al[bl . [201lt [w Kcrkwijk & Kaplan 2008[) , which 
showed that the current field strengths of the INSs 
are remarkably similar, in the range 1.0 to 3.5 x 
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10^'^ G, and that their characteristic ages of several Myr 
are substantially in excess of true ages of ~ . 5 Myr 
inferred from coo li ng and kine matics (IWalteri 120011: 
Motch et al.l l200"5l: iKaplan. van Kerkwii k. fc AndersonI 



2007t Ivan Kerkw iik fc KaplanI I200S: .Motch et al.l 120091: 
Tetzlaff et al.l l2010. 20 lH). The long periods and char- 



acteristic ages follow naturally if the INSs initially had 
much stronger fields and thus faster spindown, and 
the similar current field strengths can be understood if 
fields stop decaying at a co mmon value. Indeed, this 
was pr edicted theoretically bv lPons. Miralles. fc Gepper^ 
(|2009l ): for initially weak magnetic fields, field de- 
cay leads to only a f actor ^ 2 change that is es- 
sentially unnoticeable ('P opov et al.l l2010t ). while for 
fields above 10^'^ G field decay becomes increasingly 
important, with predicted final fields that are al- 
ways a few X 10^"^ G, independent of the initial val- 
ues. The field-decay induced heating may also help ex- 
plain the observed preponderan ce of INSs compared to 
"normal" middle-aged pulsars (iKaplaii fc van KerkwiikI 
[2009i [Faiicher-Giguere fc Kaspill2006l: IKaspill2010fr 

If the above is correct, the progenitors of the INSs 
would have been neutron stars with fields ranging from 

5 to 50 X 10^'^ G. The most magnetized would corre- 
spond to magnetars, but many would not be so ener- 
getic. Still, as their fields decay on ~ 100 kyr timescales, 
they should be brighter than expected from simple cool- 
ing, and there should be a population of relatively long- 
period, strong-field sources that are anomalously hot. 
Verifying this might not only confirm the hypothesis for 
the INS population, but also yield clues to their emis- 
sion: with the higher temperatures and stronger mag- 
netic fields, one might expect the surface to have a dif- 
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Fig. 1. — P-P diagram, with known pulsars shown by dots, and 
strong-field populations highlighted (isolated neutron stars: red 
circles; magnetars: green diamonds; strong-field pulsars with X- 
ray observations: magenta stars [filled: detected; open: upper limit 
only]). Sources referred to in the text are labeled, as are lines of 
constant magnetic field, lines of constant characteristic age, and 
an approximate death line. Overlaid on the diagram are the pre- 
dictcd spin-down tracks from the field decay model of pons et al.; 
(2009) (solid lines: dipole only; dashed lines: including a toroidal 
component with a strength of 50 times that of the dipole). The 
initial dipole fields are lolS.lS-S.w.W-S.lS q g^j- j-fjg p^jg^ ^jjgj.g 
are cross-ticks at (0.01, 0.03, 0.1, 0.3, 1) Myr. We assumed initial 
periods of 0.1 s, and divided the field at the pole by two to get the 
field at the equator for comparison with the spin-down field esti- 
mates. For low field strengths, field decay is not important, and 
hence the tracks essentially follow constant-field lines (shown by 
the diagonal thin solid lines), and at any moment the character- 
istic age (diagonal thin dotted lines) is a reasonable estimate for 
the true age. For strong fields, however, field decay is important, 
and the field and characteristic age inferred at any moment are not 
representative for the initial field and true age. 

ferent state and the spectra to show different spectral 
features that could be contrasted to the INSs. 

Checking the P — P diagram for possible INS pro- 
genitors (Fig. [ij, one finds four radio pulsars, one ro- 
tating radio transient (RRAT), and one relatively low- 
field anomalous X-ray pulsar (AXP); all have periods 
longer than 3 s and fields in excess of 3 x 10^"^ G. Of 
course, the AXP (IE 2259-1-586), like other magnetars, 
is thought to be powered by magnetic field de cay; it cer- 
tainly is anomalously hot, with kT ~ 0.4 keV (jZhu et al.l 
I2008l i. Of the others, the RRAT and one of the pul- 
sars have X-ray counterparts as well, with thermal spec- 
tra and infe rred t emperatures of 140 — 190 eV (PSR 
J1718-3817,|K^m & McLaughlin 2005; Zhu et al.ll20TTI : 
PSR J1819-1458. rMcLaughlin et al.„2007l) . These tem- 
peratures are, as the authors mention, well above expec- 
tations from simple cooling (and also well above th e limit 
of~70eV inferred for PSR B0154+61 ( Gonzalez_eLalJ 
[200l . a young pulsar with a somewhat weaker in- 
ferred field; see Fig. [T]). The best-studied source, 
PSR J1819— 1458, also has a clear absorption feature 
at 1 keV, well above the energies at which absorption is 
seen in the INSs. Among the three remaining sources, 
PSR J1814-1744 and PSR J1847-0130 have unpub- 
lished, reasonably long XMM observations. We checked 
these and found no counterparts, but as both pulsars are 
distant and extincted, this is not unexpected. Indeed all 
of these objects are considerably further than the INSs, 
and while part of this can be explained by their small 
ages leading to lower space densities, part comes from 



the narrowly-beamed/intermittent radio emission that is 
used to detect them (compared to omni-directional soft 
X-rays). 

The last source, PSR J0726-2612 (hereafter PSR 
J0726) is the subject of this paper. This 3.44 s pul- 
sar was discove red in the course of the Parkes High- 
Latitude survey (|Burgav et al.ll2006D but surprisingly has 
not seen any X-ray follow-up, despite its inferred mag- 
netic field of -B = 3 X 10^^ G and characteristic age of 
only r = P/2P — 200 kyr. Yet, it is arguably the most 
interesting, since it should be the least extincted — both 
because it has the lowest dispersion measure (DM) of all 
(70 vs. 200 to 800 cm~'^ pc for the other four pulsars) and 
because it is at what is, for a young pulsar, high Galactic 
latitude (|6"| = 4.7° vs. < 0.22° for the others). Fur- 
thermore, it may well be the closest. From its dispersion 
measure, co mbined with a model of the Galactic electron 
distribution (jCordes fc Laziol [20021 ) . one infers a rough 
distance of 3 kpc. But at that distance, its height above 
the Galactic plane is 250 pc, which is somewhat improb- 
able. A more likely alternative is that i t is at ^ 1 kp c, sit- 
uated (and born) in the Gould Belt ( Popov et al.l i2005l. 
which this line of sight passes through (and which may 
influence the dispersion measure). Below, we will scale 
our distances to dkpc = d/lkpc. 

Overall, PSR J0726 seems a prime candidate for 
comparison with the INSs, being relatively close and 
unabsorbed, and apparently intermediate between the 
stronger-field pulsars that look most like INS progen- 
itors and the weak-field normal pulsars that form the 
bulk of the population. Intriguingly, we found a pos- 
sible counterpart to PSR J0726 in the RASS: IRXS 
J072559.8-261229, with 0.027 ± 0.011 ct s^^. The nomi- 
nal separation of 112" from the pulsar is relatively large, 
but in the RASS image, the source seemed somewhat 
extended, encompassing the pulsar. Here, we present a 
Chandra X-ray Observatory of IRXS J072559.8-261229 
that confirms the identification with PSR J0726, and al- 
lows a first comparison between PSR J0726 and the INSs, 
opening the way for detailed foUowup. 

2. OBSERVATIONS & ANALYSIS 

We observed PSR J0726 with the Advanced CCD 
Imaging Spectrometer (ACIS; iGarmire et al.l 120031 ) 
aboard the Chandr a X-ray Observatory {CXO; 
'Weisskopf et alj [2000f) on 2011 June 15 for 17.9 ks 
(ObsID 12558). For best sensitivity at low energies, 
we used the back-illuminated S3 CCD, selecting the 
1/8-subarray read-out mode (with 0.4 s samphng) 
to resolve the 3-s pulse period. We processed the 
level- 1 event lists following standard procedures (using 
chandrajrepro from CIAO v4.3 and CALDB v4.4.5). 
The data show one clear source, with a J2000 position 
(from celldetect) a = 07'^26'"08?14, 5 = -26°12'38'.'7 
with an uncertainty dominated by the boresight error 
of ~ 0'.'6 (90% confidence). This is consistent with the 
radio position of PSR J0 726 in the ATNF pulsar catalog 
(|Manchester et alJ a = 07h26™08!12 ± 0!04, 

5 = -26°12'38'.'1±0'.'8. 

For our further analysis, we selected 1179 source events 
from a circular region with a radius of 5 pixels (2'.'46), en- 
compassing > 95% of the source photons, and 326 back- 
ground events from the rest of the detector area (a fac- 
tor of 1200 larger in area), with energies between 220 eV 
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Fig. 2. — Lower panel: X-ray (blue) and radio (red) pulse profiles 
of PSR J0726— 2612. The X-ray data are from our new Chandra 
observations, while the radio data are from the ATNF pulsar data 
archive with an arbitrary flux offset and scaling (these are from a 
152-s observation centered at 1374 MHz with 288 MHz bandwidth, 
taken on 2005 March 26). Both pulses are repeated twice for clarity. 
The uncertainty on the X-ray TOA is 0.013 cycles, but in compar- 
ing that to the radio TOA we must include the uncertainties on 
the ephemeris which contribute substantially. The apparent shift 
between the radio and X-ray pulses is 0.07 ± 0.11 cycles, as shown 
by the black error bars. Upper panel: the median energy in each 
phase bin. 

(the recommended lower-bound for data from the ACTS; 
any lower and the calibration becomes unreliable and 
soft flares become increasingly dominant) and 1.1 keV 
(where the effects of higher-energy flares are minimized; 
the source is barely detected above 1 . 1 keV) . We cor- 
rected the event times to the Solar System barycenter 
using axbary, assuming the radio position. 

2.1. Timing Analysis 

To look for X-ray pulsations, we first determined the 
frequency that ma ximized the power in a Zf periodogram 
(Rayleigh statistic; IBuccheri et al.llT983l) at different mul- 
tiples of the radio frequency. We found no obvious peak 
at the expected frequency of the pulsar {Zf = 0.7), but a 
strong one {Zf = 21.7) at roughly double the frequency, 
0.290493 ± 0.000( 302 Hz (where the uncertainty was cal- 
culated following lRansomll2001[ ). We found no evidence 
for significant power in higher harmonics {Zf = 1.2 and 
0.70 at the third and fourth harmonic, respectively). We 
folded our data on half our best-fit frequency to construct 
a binned light-curve (with 16 phase bins; Figure [2]) • As 
expected from the lack of other harmonics, a sinusoid 
provided a good fit (x^ — 11-7 for 13 degrees of free- 
dom [dof]), with a pulsed fraction (semi-amplitude) of 
27% ± 5%. The implied TOA, defined as the maximum 
of the sinusoid closest to the middle of the observation, 
is MJD 55727.6883125 ± 0.0000006. The ephemeris in 
the ATNF catalog (based on data from 2005-2006) pre- 
dicts a frequency of 0.29049679119 ± 9.6 x 10"^° Hz at 
the epoch of our observation, which differs by ^ 2 cr from 
our measurement. Further radio observations would help 
determine whether this is due to chance or due to some- 
thing intrinsic (a glitch or timing noise). 

In order to compare the phase of the radio and X- 
ray pulses, we retrieved an archived Parkes observation 
of PSR J0726 from 2005 March 26. The pulse pro- 
file, shown in Fig. [2l was substantially noisier than that 



shown in iBurgav et al.l ()200(]D , but still sufficed to de- 
termine a TOA. From the highest point in the pro- 
file, we find MJD 53455.29291874 ± 0.00000004, with 
the uncertainty being ±1 phase bins. With this TOA 
and the ATNF ephemeris, we infer a phase difference of 
TOAradio - TOAx = 0.07 ± 0.11 cycles, where the un- 
certainty is dominated by uncertainties in the ephemeris 
(see Figure [2]) . 

2.2. Spectral Analysis 

Spe ctral fitting was done in sherpa (jRefsdal et al.l 
|2009( ). For the fits, we only included events between 
0.32-1.1 keV, as the response below 320 eV is not well 
understood. This leaves 1003 source and 267 background 
counts. We follow recommendations in modeling the low- 
energy background with a power-law instead of subtract- 
ing it (although given that the expected background in 
the source area is less than 1 count, it does not influence 
the results) . We binned the source counts such that each 
bin contained at least 25 counts and was at least 29 eV 
wide, so that statistics are a good approximation, 
and we do not oversample the instrumental resolution of 
~100eV. 

We tried three main source models: a 
thermal blackbody, a non-thermal power 
law, and a simple neutron-star atm os phere 
model (NSA: IPavlov. Shibanov. fc ZavlinI 11991 
iZavlin. Pavlov, fc Shibanovl |199( ?)" all modified by 
interstellar absorption (using xsphabs; the parameters 
changed only minutely with other choices, such as 
xstbabs and xswabs). We tried NSA models with and 
without a strong magnetic field; for both, we fixed the 
mass and radius of the neutron star at 1.4 and 
10 km, respectively. Our results are summarized in 
Tabic [T] (there and below, temperatures and radii refer 
to those measured by a distant observer). 

The NSA models gave reasonable fits, but implied im- 
plausibly small distances: 33l^g and 160lgj!)°pc for the 

B = and B = 10^"^ G models, respectively. A power- 
law model gave a significantly worse fit than the ther- 
mal models (x^ = 76.9 vs. 45.1 for 67 dof). We found 
that, similar to the INSs, a simple blackbody gave a good 
fit (Fig. [3]) and the most reasonable parameters In 
what follows, we use that as our baseline model. We find 
there is significant covariance between the fitted param- 
eters, which, within la, leads to changes in absorption of 
4 X 10^*^ cm^^ and in temperature of ~5eV. To explore 
the influence of the absorption column, we also tried flts 
where we held A^h fixed at 1 or 10 times the dispersion 
measure (70cm~^pc = 2.1 x lO^^cm"^). These fits were 
somewhat worse than those with A^h free, but still had 
reduced less than 1. 

To constrain any additional non-thermal emission, we 
tried adding a power law to the blackbody. We found it 
had negligible effect on the parameters, and, assuming 
r = 3, we infer a la upper limit of the photon rate of 

^ With the blackbody model and the observed 0.32- 
l.lkeV count rate of 0.056s^^, the expected count rate 
for ROSAT PSPC is 0.027 s'^ (estimated using the 
Portable, Interactive Multi-Mission Simulator [PIMMS]; 
http: //heasarc . gsfc .nasa.gov/Tools/w3piiirais .html ). Thus, 
most likely PBK J0726 is responsible for all of the flux of IRXS 
J072559. 8-261229. 
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Fig. 3. — The best- fit blackbody model along with residuals 
(lower panel). Inset: best-fit contours of kT vs. normalization 
R/d, with 1-, 2-, and S-cr joint confidence contours (Ax^ = 
2.30, 6.1711.8). The best-fit value is the red circle, while the blue 
square and black diamond are for the combined minima and max- 
ima, respectively. 

7 X 10~^ cm^^ keV~^ at 1 keV, corresponding to < 3% 
of the flux in the 0.32-1.1 keV band. 

We also tried fitting the spectrum of PSR J0726 
with a model that included an additive Gaussian line, 
in order to constrain the presence of absorption fea- 
tures similar to those seen for the INSs , which have 
equivalent widths o f up to several 100 eV (jHaber II 120071 : 
Ivan Kerkwiik fc Ka plan 2007|). We fixed the full width 
at half maximum (FWHM) of the line to 0.1 keV, com- 
parable to the energy resolution of the back-illuminated 
ACIS-S3 CCD, and varied the line energy between 
0.3 keV and 1.1 keV, fitting only for its amplitude (allow- 
ing both positive and negative amplitudes). As expected 
from Fig. |3l we did not find evidence for any absorp- 
tion line. The largest improvement in the fit occurs near 
0.6 kcV, for an emission line with a height of 35%, but 
even this reduces by only 4, which is not significant 
given our number of trials (based on simulations, we es- 
timate a false-alarm probability of ^^30%). We infer a 
68%-confidence upper limit of 50 eV to the equivalent 
width of any feature. 

Finally, we looked for phase-resolved trends in two 
ways. First, we determined the median energy as a func- 
tion of pulse phase (upper panel in Figure [2]), finding lit- 
tle change except perhaps a hint of hardening at the max- 
imum that is coincident with the radio peak, although 
given the current statistics this is not significant (the 
for the median energy is 9 for 15 degrees-of- freedom). 
Second, we separated the data into two bins: one con- 
taining the two maxima and one the two minima, each 
with 50% of the exposure time. From our fit (Table[T]and 
Figure|3]), we find negligible difference in temperature be- 
tween the maxima and minima, consistent with the lack 
of change seen in median energy. Thus, the pulsations 
appear to reflect mostly changes in emitting area. 

2.3. Spatial Analysis 

To check for possible extended emission, e.g., from a 
pulsar wind nebula, we compared our observations with 
a simulated point-spread function (PSF), created us- 
ing ChaRT (the Chandra Ray Tracer) and MARX (version 
4.5.0), for a spectrum consistent with our best-fit black- 
body. We find good agreement for radii > 2 pixels (1"): 
although with only w 180 counts at these radii, the con- 



straints are not strong. Oddly, the two bins at the small- 
est radii, 0-0.5 pixels and 0.5-1.0 pixels, had some devi- 
ations, with the model PSF underpredicting the counts 
in the first bin by about 6 a (30%), and overpredicting 
the second by about 4cr (26%). Since it is unlikely that 
the pulsar is more centrally concentrated than a point- 
source, we probably are simply reaching the limits in the 
model (we experimented with the Chandra processing, 
turning off pixel position randomization, but found that 
this made little difference). We conclude that there is no 
evidence for extended emission on scales of > 1". 

3. DISCUSSION 

Regardless of the path by which PSR J0726 and the 
INSs got to their present states, we can examine how 
closely their current properties resemble each other and 
whether the unique information from either helps under- 
stand the thermal emission from this class of neutron 
stars as a whole. Furthermore, we can see how PSR 
J0726 might help us understand the evolution of high-_B 
pulsars and the origin of the INSs. We consider these in 
turn. 

3.1. Comparison with the INSs 

PSR J0726 was selected to be similar to the INSs in 
having a long spin-period, P > 3 s, and strong magnetic 
field, P = (1 - 3) X 10" G. With a thermal spectrum 
with kT = 87 ± 5 eV, and i? « 6 km (at the nominal 
distance of ~ 1 kpc in the Gould Belt), its X-ray prop- 
erties are similar as well (indeed, it would be classified 
as an INS based on these). Its spectrum has no evi- 
dence for absorption, u nlike RX J1308. 6-1- 2127, which 
has strong absorption (jHaberl et al.l I2003D while hav- 
ing a very similar magnetic field and a slightly higher 
tempe rature {B = 3.4 x 10^'^ G, ^apla n & van Kerkwiik 
"2005^; kT = 102 eV. ISch^ope et al.M2007D . However, 
RX J0720.4-3125, another INS with a similar ff slightly 
lower field {B = 2.4 x 10" G, iKaplan fc van Kerkwii'kl 
l2005aD . had weak absorption, with an equivalent width 
of ~ 5eV, which we would not be able to detect, 
when its temperat ure was most si milar to PSR J0726 
{kT = 86 eV, .Haberl et al.l I2004D . Intriguingly, this 
source changed, developing much stronger absorption 
while becoming hotter (kT = 94eV.lde Vries et alll2004l: 
iHaberl et all 120061: Ivan Kerkwiik et alj \206W . Taking 
our inferred temperature at face value, it is thus pos- 
sible that PSR J0726 emits identically to an INS of the 
same temperature and dipole magnetic field strength. If 
so, at higher sensitivity it may show weak absorption at 
~ 300 eV. (We caution, however, that the dependencies 
of absorption energy and strength on temperature and 
magnet ic field strength remain puzzlin g for the INSs as 
a class; iKaplan &: van Kerkwiikll2009al lbl.'l 

The inferred luminosity Lx = 1-5 x lO'^^d^p^. ergs~^ 
(0.32-1.1 keV) is about half the rotational energy loss 
rate E = Att^IP/P^ = 2.8 x lO^^ergs^^ (for a moment 
of inertia / — lO'^^gcm^). This is on the low side of 

what is seen for the INSs (Lx ~ (0.5 . . . 100)£'), possibly 
just a consequence of PSR J0726 having a shorter period 
for the same B (for dipole spin-down, E oc B'^/P^). It 
is much larger than the non-thermal emission typically 
observed for ro tation-powered pulsars {Lx ~ 10"'^ -E, 
iBecker fc Triimper. .1997.') . although we cannot exclude 
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TABLE 1 

Results of X-ray fits to PSR J0726-2612 



Model 


Nn 


r or kT 


A or R/d'' 










(10^0 


(N/A 
or eV) 


(lO^-'^s-icm-^keV"^ 
or kmkpc"^) 


(crgs^-"- cm^^) 


(erg s^^ cm^^) 





Blackbody 


8+4 


87 ±5 




4.5 


12.9 


45.1/67 


Blackbody 


2 


95 ± 2 


Ci.O_Q 2 


4.5 


5.6 


48.1/68 


Blackbody 


21 


74 ±2 


i7.it}i 


3.6 


33.8 


55.0/68 


NSA {B - 0) 


141^ 


22 ±3 


3071?™ 


4.2 


22.7 


52.8/67 


NSA (B = 10" G) 




42 ±4 




4.3 


19.2 


52.8/67 


Power Law (P free) 


2011 


6.4 ±0.4 


7 6+1° 

'•"-0.8 


4.2 


60.6 


76.9/67 


Power Law (P = 3) 


0+0.1 


3 


9-OlSi 


3.2 


3.2 


195/66 


Blackbody phase-resolved: 


9tl 










Minima'^ 


85 ±6 


5.81?;? 


3.3 


8.5 


24.3/37 


Maxima'^ 


9tl 


86 ±6 


7 1+3.5 
'■^-2.2 


5.1 


13.2 


47.1/51 



Note. — Quantities without uncertainties were held fixed during the fit. All uncertainties are for l-cr confidence, and all other parameters 
were allowed to vary during the calculation of the uncertainties. 

For the NSA models, the mass is fixed to 1.4 Mq and the radius is fixed at 10 km, so the uncertainties calculated for R/d are based on 
the best-fit distance. The normalization of the power-law model A is the photon rate at an energy of 1 keV. 
^ The fluxes are given in the 0.32-1.1 keV band and are given both as observed and corrected for absorption. 

Phase-resolved spectroscopy using a blackbody model. The two phase bins each include 50% of the data, and were flt simultaneously 
with the same absorption. 



that PSR J0726 has a shnilar non-thermal component: 
our Icr upper hmit of 3% of the 0.32-1.1 keV flux for the 
contribution of a power-law with F = 3 (a value typical 
for pulsars) implies Lpl < O.Olbd^^^E. 

The pulse properties of PSR J0726 also echo those of 
the INSs, which show smooth profiles dominated by the 
fundamental or the first harmonic. The pulse fraction of 
27 lb 5 % is a little la rger than the ^1 to 18% seen for the 
INSs (jHaberll 120071) ■ but it is not dramatically different. 
We note that without the radio period, we might have 
actually identified the pulse period of PSR J0726 as half 
of the true value. The lack of temperature variation is 
somewhat strange (cf. up to 10% changes in kT over the 
pulse of RX J1308. 6-1-2127, for example; iSchwope et al.l 
I2005D . as, naively, it implies that we are seeing changes 
in projected area of large regions of similar temperature 
with anything in between so much colder that it does 
not contribute. It may be that instead the whole surface 
emits at more or less a uniform temperature, and that the 
pulsations reflect anisotropics in the emission, perhaps 
associated with the strong magnetic field. The same may 
hold for the INSs, which typically also show only modest 
changes in temperature with pulse phase. 

A possible difference with the INSs is that PSR J0726 
has an order-of-magnitude smaller characteristic age. 
This may again be a consequence of its shorter period 
for the same B (for dipole spin-down, r oc (P/B)'^), but 
might also point to an evolutionary difference (see be- 
low). Another characteristic of INSs is that they have 
faint optical counterparts, but with optical/UV fluxes 
a factor of 6-50 above the extrapolation of the X-rays. 
This is in contrast to middle-aged pulsars such as PSR 
B0656-I-14 and Geminga, which have ultraviolet fluxes 
that are more consistent with the X-ray ext rapolations 
(|Shibanov et al.ll2005l : iKargaltsev et al.ll2005l ). It will be 
interesting to see whether PSR J0726 has an excess or 
not. 

Of course, the most glaring difference between PSR 
J0726 and the INSs is the radio emission. Numerous 



searches have yet to find confir med radio emission (co- 
herent or bursty) from the INSs (|Kondratiev et al.ll2009L 
and references therein). This may be a result of their 
location near the pulsar "death line," or more simply 
a consequence of their narrow radio beams: the 50% 
width of the radio pulse from PSR J0726 is < 1% of 
the pulse period (consistent with general trends for long- 
period pulsars) so the chance of missing the radio pulse 
is large. 

Overall, we conclude that with the exceptions of its ra- 
dio emission and characteristic age, the properties of PSR 
J0726 agree with those of the INSs. Based o n the com- 
pilatio ns of iKaplan fc: van Kerkwiilj (f2009bD . iZhu et ahl 
(|2011| ) , and references therein, it is the best pulsar analog 
to the INSs to date. 

3.2. Implications for the Evolution of High-B Neutron 

Stars 

In order to look for evolutionary trends, we can com- 
pare PSR J0726 to sources with similar properties in the 
P-P diagr am, as wel l as with the evolution expected from 
the work of fPons et al. (2009), which seems to explain the 
origins of the INSs (Fig. [1]). A first comparison is with 
INSs of similar field and temperature, RX J0720.4-3125 
and RX J1308. 6-1-2127. As mentioned, the main dif- 
ference is that the characteristic age of PSR J0726 is 
200 kyr, while those of the INSs are > 2 Myr. This could 
imply that the true age of PSR J0726 is closer to its char- 
acteristic age than is the case for the INSs, and maybe 
that it has underg one less i? decay. This is consistent 
with the models of iPons et all ()2009[) : 3 x 10^^ G is near 
the dividing line where field-decay becomes important, 
and the true age is expected to be no more than a factor 2 
shorter than the characteristic age (see Fig. [ij. Taking 
the models at face value, one infers that PSR J0726 was 
born with a dipole field of just below 10^"* G, while most 
INSs would have had a field about a factor 2 stronger. 

As a second comparison, we can compare PSR J0726 
with two X-ray-detected objects that have almost the 
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same period: the modest-g INS RX J04 20.0-5022 
(1.0 X 10^^ G, Kaplan fc van Kerkwiid 120111 ) and the 
high-B pulsar PSR J1718-3817 (7.4 x lO^^ |Zhu et al.l 
1201 ID . The former is very similar to PSR J0726, with 
the main difference being that it is cooler {kT w 45 eV). 
From the evolutionary tracks, it seems most likely RX 
J0420. 0—5022 initially simply had even lower magnetic 
field strength, which decayed by a factor ^ 2 or so, and 
that, like PSR J0726, its true age is at most a factor 2 
shorter than its characteristic age of 2 Myr. 

From the evolutionary tracks, PSR J1718— 3817 would 
be a good candidate progenitor for INSs like RX 
J0720.4— 3125, and indeed it has a largely thermal spec- 
trum consistent with a blackbody with a significantly 
hotter temperature {kT « 190 eV), suggestive of on- 
going field decay. However, the emitting radius of 2 km 
and pulsed fraction of ^ 50% differ from those found for 
the INSs (and PSR J0726), and are more reminiscent 
of even younger high -i? pulsars like PSR J1119— 6127 
([Gonzalez et al.ir2005f) , which also have high pulsed frac- 
tions and smaller radii (although distance and fitting un- 
certainties make the radii for all these objects poorly con- 
strained). It may be that this reflects some combination 
of ongoing B decay and non-thermal heating; all these 
younger, strong-field objects have substantially larger 
spin-down energy losses. Alternatively, the stronger field 
may lead to more anisotropic emission. 

Overall, it seems that the INSs and high-B pulsars 
form a single family (also see Kaspi"2010' for further dis- 
cussion) , with initial fields that were stronger than their 
current ones, and that caused rapid initial spin-down. 
The long periods for the slower spinning objects also 
imply low rotational energy losses and thus weak non- 
thermal emission. Combined with some additional heat- 
ing due to field decay, this makes the thermal emission 
stand out more than it would in normal pulsars. 

A continuing puzzle, however, is the difference with 
the magnetars (|Ng fc Kaspill201in . which occupy similar 
parts of the P-P diagram, but are clearly much more 
strongly affected by magnetic field decay. It may be 
that those were born with much stronger toroidal com- 
ponents to the field (which affect the spin-down only in- 
directly, via the more rapid field decay). In any case, 
there may be a smooth continuum, given objects like 
PSR J1846-0258 {B = 4.9 x 10" G) wh ich exhibited 
a sudden, magnetar-like X-ray outburst (iGayriil et al.l 
[2008tlKumar fc Safi-Harbl[2008l: iNget al.ll2008n . the INS 
RX J0720.4— 3125 that exhibited a possible magnetic 
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recon figuration (|de Vries et al.|[200^ Ivan Kerkwiik et all 
I2007f) . the transient mag netar XTE J1810-1 97 that 
was detected in outburst ijibrahim et al.l I2004D but in 
quies cence appears more like an INS (jBernardini et al.l 
120111 ). or the putative low-B magnetar SGR 0418- f 5729 
(< 8 X 10^2 Q. IRea et al.ll2010l : iTuroUa et alllMl . 

4. OUTLOOK 

Our Chandra observations show that the relatively 
nearby high-i? pulsar PSR J0726— 2612 has properties 
similar to those of the INSs, showing analogous ther- 
mal emission and similar smooth pulsations. Because of 
its radio emission, PSR J0726— 2612 may help us clarify 
the nature of the emission and the pulsations. First, for 
many radio pulsars the observed polarization position an- 
gle curve is well described by th e rotating vector model 
(|Radhakrishnan fc Cookd 119690 , in which the position 
angle is assumed to be aligned with a dipolar magnetic 
field. Using this model, one can infer the (mis) alignment 
between the spin and magnetic axes, as well as the an- 
gle between the spin axis and line-of-sight. Using such a 
model for PSR J0726 and combining it with fits to the X- 
ray light-curve will help resolve the a mbiguities that lim- 
ited lightcurve analyses for IN Ss (Bra ie fc Rom ani 200^; 
IZane fc Turo"nall2006l : Illoll2007i) . Second, unlike for the 
INSs where only X-ray and optical astrometry (both of 
limited precision) are possible, for PSR J0726 traditional 
radio techniques are available (proper motions with the 
Very Large Array, and potentially parallaxes with the 
Very Long Baseline Array) . This means that even though 
the likely distance is ~ 1 kpc, which would be well outside 
of reach for an optical parallax, we can hope for a geo- 
metric distance measurement for PSR J0726 that would 
constrain its emitting radius. 
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